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Abstract 

We recently completed a calculation of the process e^e^ QQ + X, where 
Q is a heavy quark, at order 0(0^). As a first application of this calculation we 
compute the momentum correlations of bb pairs at next-to-leading-order. This 
quantity is interesting since it may affect the determination of Rb as measured 
in decays. We find that the next-to- leading corrections are of moderate 
size, thus confirming the conclusions that can be drawn from the leading-order 
calculation. 
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1 Introduction 

Radiative corrections to jet production in e^e~ annihilation have been known for 
a long time 0, Previous calculations were, however, performed for massless 
quarks. In most practical applications this is sufficient, since at relatively low energy, 
the b fraction is strongly suppressed, and at high energy (i.e. on the peek and 
beyond) mass effects are presumably suppressed. Nevertheless there are several rea- 
sons why a next-to-leading-order calculation is desirable. First of all, at sufficiently 
high energies, top pairs will be produced, and mass effects there are very likely to 
be important. A second important reason is to understand the relevance of mass 
corrections due to bottom production to the determination of as from event shape 
variables. As a third point, quantities such as the heavy flavour momentum correla- 
tion [Q, although well defined in the massless limit, cannot be computed using the 
massless results of refs. fl], ||, ^ . 

We have recently completed a next-to-leading-order calculation of the heavy- 
flavour production cross section in e^e~ collisions, including quark mass effects. As 
an illustration, in fig. |l| we show the Feynman diagrams for a Born term (a), a virtual 
correction term (b), and two real next-to- leading contributions (c,d). Next-to- leading 
corrections arise from the interference of the virtual graphs with the Born graphs, 
and from the square of the real graphs. Observe that we always deal with the cross 
section for the production of the heavy quark pair plus the emission of at least one 
extra particle (i.e. a gluon or a quark). The inclusion of virtual graphs with only a 
QQ pair in the final state is not needed if one only computes 3-jet related quantities. 
Furthermore, by using available results on the total heavy fiavour cross section at 
order (see and references therein), one can avoid altogether to compute this 
kind of virtual graphs. 

Virtual graphs, besides the usual ultraviolet divergences, which are removed by 
renormalization, have also infrared and coUinear divergences. These cancel when suit- 
able infrared safe final-state variables are considered. Our treatment of the infrared 
cancellation is such that the final result is expressed as a partonic event generator, in 
which pairs of weighted correlated events are produced. Any shape variables can be 
computed for each generated event. Only infrared-safe variables give rise to a finite 
distribution. No arbitrary cutoff is needed in this calculation in order to implement 
the cancellation of virtual and real infrared divergences. Therefore, one does not 
have to worry about taking the limit for a vanishing soft cutoff. Thus, this method is 
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Figure 1: Some of the diagrams contributing to the process Z^QQ + X : 
a Born graph (a), a virtual graph (h), a real emission graph (c), and a real 
emission graph with light quarks in the final state (d). 



similar to the one of ref. 0, which was used there to compute a large class of shape 
variable distributions for the LEP experiments. 

We have verified that our calculation has the correct massless limit, by comparing 
our result with that of the massless calculation of ref. for all the shape variables 
listed there, and for the moments of the energy-energy correlation, proposed as a 
benchmark in ref. 0. 

Very recently, two calculations have appeared in the literature that address the 
same problem |^ ^. They both use a slicing method in order to deal with infrared 
divergences. We were able to perform a partial comparison of our result with that of 
ref. 0, and found satisfactory agreement. In the older work of ref. |jlO|, a calculation 
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of the process e^e — * QQgg has been given, but virtual corrections to the process 



e e -^QQg were not included. In ref. |11], the NLO corrections to the produc 



tion of a heavy quark pair plus a photon are given, including both real and virtual 
contributions. 

We will give full details of our calculation in a future publication. In the present 
letter we would like to present a first application, which could not be dealt with using 
the massless results. This is the average momentum correlation r defined as 

_ {X1X2) - {xf 

{xf ' ^ ' 

where Xi(2) are the Feynman x of the produced B{B) mesons, Xi(2) = Pb{b)/p^b'^^\ 
(x) = (a;i(2)). In ref. 0] r is computed at leading order in the strong coupling constant. 
The interest in this quantity lies in the fact that it may affect the determination of 



{{1 - x,){l - X,)) = ^b+[^yc + 0{al) 



Rb when the tagging efficiency for B mesons is determined by looking at both 
the single-tag and the double-tag events. A simplified description of how this works 
is given in ref. It is convenient to rewrite r in the following way 

^_ {(l-x^)il~x,))-{l-xf _ ^^ 2) 

{xf 

We notice that, thanks to the above form, the contribution of virtual graphs with 
only a QQ pair in the final state is not needed for this calculation. Defining now 

2 

(l-x) = ^a + 0{al), (1.3) 

where as is the MS , 5-fiavour coupling constant evaluated at a scale /i equal to the 
annihilation energy Q. We have 

r = ^ & + (^) ' (c + 2a6 - a') + 0{al) . (1.4) 

In the particular case of the computation of the average momentum correlation, 
we have treated as massive the heavy quark directly produced from the electroweak 
vertex. Secondary heavy quarks produced by gluon splitting have been treated as 
light. Similarly, we have not included graphs in which a light fiavour is produced by 
the electroweak vertex, and the heavy quark pair is produced by a gluon. Although 
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the inclusion of these contributions does not pose any particular problem, we believe 
that this quantity is the closest to what is usually defined to be the correlation induced 
by hard radiation. Other dynamical correlation mechanisms, arising for example from 
the production of four heavy flavours, or from the indirect (e.g. via gluon splitting) 
production of a heavy flavour pair in an event with primary production of light quarks 
can be examined separately. Their influence upon the determination of Ti, depends 
upon a more detailed specification of the experimental setup, such as, for example, 
on whether cuts are applied that remove soft heavy quarks. 

Results for r are displayed in table |l[ From the table it is apparent that r cannot 



m 


a 


b 


c 


r 


1 GeV 


12.79(1) 


0.6628(1) 


155.5(3) 


0.1055a, + 0.23(1)^2 


5 GeV 


7.170(2) 


0.6182(1) 


50.94(5) 


0.0984a, + 0.213(1)^2 


10 GeV 


4.858(1) 


0.5432(1) 


26.23(2) 


0.0865 a, + 0.2004(6)a2 



Table 1: Mass dependence of the coefficients a, 6, c and r. The digit in parenthesis 
(to be taken as zero if not given) represents the accuracy of the last digit. 



be computed in the massless limit, since the quantities c and a are both plagued by 
coUinear divergences. The coefficients of the expansion of r itself do instead converge 
in this limit. This fact was discussed at length in ref. [Q. From the table we also see 
that the radiative corrections to r are small. Assuming as{Mz) = 0.118 (correspond- 
ing to the PDG average |]13[), we have in leading order r = 0.0984 x a, = 0.0116, 
and in next-to-leading order r = 0.0984 x a, + 0.213 x a^ = 0.0146. Thus, radiative 
corrections are reasonably under control, and do not spoil the main conclusion of 
the leading-order calculation. With the assumption of a rough geometric growth of 
the expansion, we can give an estimate of the theoretical error due to higher orders: 
r = 0.0146 ± 0.0007. We have also computed the quantity r', defined in ref. Q as 

, (a;ia;2cut(l,2)) - {xf 

r = ^ (1.5) 

{xf 

where cut(l,2) is defined to be 1 if the quark and antiquark are in opposite hemi- 
spheres with respect to the thrust axis, and zero otherwise. We define 

fy / Ol 

{1-XI-X2 + xiX2 cut(l, 2)) = 6' + c' + 0{al) , (1.6) 

27r \2tt/ 
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and obtain 



(c' + 2ah' - a^) + 0{al) . (1.7) 
The quantities 6', d and r' are given in table Assuming as before as{Mz) = 0.118 



b'+ , , 
27r \27tJ 



m 


h' 


d 


d 


1 GeV 


0.3800(7) 


158.4(3) 


0.0605(l)a, + 0.12(1)^2 


5 GeV 


0.3653(6) 


50.9(1) 


0.0581(l)a, + 0.119(2)^2 


10 GeV 


0.3423(5) 


25.51(6) 


0.0545(l)a, + 0.133(1)^2 



Table 2: Mass dependence of the coefficients a, 6', d and r' . 



we get r' = 0.0064 at leading order, and d = 0.0083 at next-to-leading order. 

As far as its perturbative expansion in powers of is concerned, the average 
momentum correlation is a quantity that is well behaved in perturbation theory, 
and it is also quite small. Since its effect is typically of the order of 1%, one may 
worry that non-perturbative effects, of order A/Q (where A is a typical hadronic 
scale) may compete with the perturbative result. This is a very delicate problem, 
since we know very little about the hadronization mechanism in QCD. In ref. [Q, 



this problem was addressed in the context of the renormalon approach to power 
corrections. It was shown there that (in the renormalon approach) corrections to 
the momentum correlations are at least of order (A/Q)"^, and thus negligible at LEP 
energies. Although this result cannot be considered as a definitive answer to the 
problem, it is at least an indication that power corrections to this quantity are small. 



References 

[1] R.K. Ellis, D.A. Ross and A.E. Terrano, Nucl. Phys. B178(1981)421. 

[2] J. A.M. Vermaseren, K.J.F. Gaemers and S.J. Oldham, Nucl. Phys. 
B187(1981)301. 

[3] K. Fabricius, G. Kramer, G. Schierholz and I. Schmitt, Z. Phys. Cll(1981)315. 

[4] R Nason and C. Oleari, Phys. Lett. B387(1996)623, hep-ph/9607347. 

^Small numerical discrepancies with ref. [|| are due to round off errors, and to the inclusion of 
the electromagnetic term in the present formulae. 



-7- 



[5] K.G. Chetyrkin, R. Harlander, J.H. Kiihn and M. Steinhauser, preprint MPI- 
PHT-97-012, hep-ph/9704222. 

[6] Z. Kunszt and P. Nason, "QCD", in "Z Physics at LEP 1", eds. G. Altarelli, 
R. Kleiss and C. Verzagnassi (Report CERN 89-08, Geneva, 1989). 

[7] P. Nason and B.W. Webber, "QCD", in "Physics at LEP2", eds. G. AltareUi, 
T. Sjostrand and F. Zwirner (Report CERN96-01, Geneva 1996). 

[8] G. Rodrigo, Nud. Phys. Proc. Suppl. 54A(1997)60, hep-ph/9609213; 
G. Rodrigo, Ph. D. Thesis, Univ. of Valencia, 1996, hep-ph/9703359; 
G. Rodrigo, A. Santamaria and M. Bilenkii, preprint FTUV-96-80, Mar 1997, 
hep-ph/9703358. 

[9] W. Bernreuther, A. Brandenburg and P. Uwer, preprint PITHA-97-5, Mar 1997, 
hep-ph/9703305. 

[10] A. Ballcstrero, E. Maina and S. Moretti, Phys. Lett. 294(1992)425; Nud. Phys. 
415(1994)265. 

[11] L. Magnea and E. Maina, Phys. Lett. B385(1996)395, hep-ph/9604385. 

[12] ALEPH CoUab., R. Barate et ai, preprint CERN-PPE-97-018, Feb 1997; 
ALEPH CoUab., R. Barate et al, preprint CERN-PPE-97-017, Feb 1997; 
ALEPH CoUab., D. Buskuhc et aL, Phys. Lett. B313(1993)535; 
ALEPH CoUab., D. Buskuhc et aL, Phys. Lett. B313(1993)549; 
DELPHI CoUab., P. Abreu et al, Z. Phys. C70(1996)531; 
DELPHI CoUab., P. Abreu et al, Z. Phys. C66(1995)323; 
DELPHI CoUab., P. Abreu et al, Z. Phys. C65(1995)555; 
L3 CoUab., O. Adriani et al, Phys. Lett. B307(1993)237; 
OPAL CoUab., K. Ackerstaff et al, Z. Phys. C74(1997)l; 
OPAL CoUab., R. Akers et al, Z. Phys. C65(1995)17; 
SLD CoUab., K. Abe et aL, Phys. Rev. D53(1996)1023. 

[13] Particle Data Group, R.M. Barnett et al., Phys. Rev. D54(1996)l. 

[14] P. Nason and B.R. Webber, Phys. Lett. B395(1997)355, hep-ph/9612353. 



